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ABSTRACT 

We investigate the properties of "star forming regions" in a previously pub- 
lished numerical simulation of molecular cloud formation out of compressive mo- 
tions in the warm neutral atomic interstellar medium, neglecting magnetic fields 
and stellar feedback. We study the properties (density, total gas-|-stars mass, 
stellar mass, velocity dispersion, and star formation rate) of the cloud hosting 
the first local, isolated "star formation" event and compare them with those of 
the cloud formed by the central, global collapse event. In this simulation, the 
velocity dispersions at all scales are caused primarily by infall motions rather 
than by random turbulence. We suggest that the small-scale isolated collapses 
may be representative of low- to intermediate-mass star-forming regions, with gas 
masses (Mgas) of hundreds of solar masses, velocity dispersions (T„ ~ 0.7 km s~^, 
and star formation rates (SFRs) ~ 3 x lO^^Moyr"^, while the large-scale, mas- 
sive ones may be representative of massive star forming regions, with Mgas of 
thousands of solar masses, (T„ ~ a few km s~^, and SFRs ~ 3 x lO^^Mgyr"^. 
We also compare the statistical distributions of the physical properties of the 
dense cores appearing in the central region of massive collapse with those from 
a recent survey of the massive star forming region in the Cygnus X molecular 
cloud, finding that the observed and simulated distributions are in general very 
similar. 

However, we find that the star formation efficiency per free-fall time (SFEg) of 
the high mass region, similarly to that of OMC-1, is low, ~ 0.04. In the simulated 
cloud, this is not a consequence of a "slow" SFR in a nearly hydrostatic cloud 
supported by turbulence, but rather of the region accreting mass at a high rate. 
Thus, we find that measuring a low SFEg may be incorrectly interpreted as 
implying a lifetime much longer than the core's local free-fall time, and an SFR 
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much slower than that given by the free-fall rate, if the accretion is not accounted 
for. We suggest that, rather than requiring a low value of the SFEg everywhere 
in the Galaxy, attaining a globally low specific SFR requires star formation to 
be a spatially intermittent process, so that most of the mass in a GMC is not 
participating of the SF process at any given time. Locally, the specific SFR of a 
star-forming region can be much larger than the global GMC's average. 

Subject headings: ISM: clouds — ISM: evolution — Stars: formation 
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Introduction 



The formation of m assive stars is currently a matter of intense debate (e.g. 



Zinnecker &: Yorke 



20071 ). High-mass star forming regions are characterized by more 
extreme physical conditions than their low-mass counterparts, containing cores of size, 
mass, and velocity dis persion roughly an order of raagnitude larger than those of cores in 



low-mass regions (e.g. 



1999 



Kurtz et al. 



Jiiina. Mvers fc Adams 



2000 



Beuther et al 



1999 



Lee fc Myers 



1999 



Garay fc Lizano 



20071 ). In particular, typical values of the properties 



of clumps within high-mass star forming regions (HMRs) are sizes 0.2-0.5 pc, mean densities 
n ~ 10^ cm~^, masses between 100 and 1000 M©, and velocity dispersions ranging between 
1.5 and 4 km s~^. In turn, the clumps break down into even denser "cores" that are 
believed to be the immediate precursors of single or gravitationally bound multiple massive 
protostars. 

The high velocity dispersions of clumps in HMRs are generally interpreted as strong 



turbulence that man ages to support the clumps against gravity (e.g.. 



McKee fc Tan 



Garay &: Lizano 



1999 



20031 ). However, the notion of "turbulent support" is difficult to maintain at 



the scales of these cores. If turbulence is to "support" a cloud against self-gravity, it must 
consist of isotropic motions, with typical scales significantly smaller than the size of the 
cloud, so that collectively thes e motions can act as a source of isotropic pressure, similar ly 



to thermal molecular motions (IMac Low fc Klessen 



2004 



Ballesteros-Paredes et al 



2003). 



However, turbulence is a fiow regime i n which the largest velocity differences are 



associated with the largest separations (e.g. 



Frisch 



19951 ). so that, in a turbulence- induced 



density fiuctuation (i.e., a cloud or clump), the largest velocity difference is expected to 
occur at scales comparable to the size of the clump itself, thus being highly anisotropic with 
respect to it. This property is indeed observed in molecular clouds, in which a principal 
component analysis of the velocity structure shows that the dominant mode is dipole-like 
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( iHeyer &: BruntI 120071 ). indicating either shear or compression, but not sohd-body rotation 
(M. Heyer, 2009, private communication). Being dipolar, the principal velocity component 
appears unsuited to act as a supporting form of pressure capable of opposing self-gravity. 



Moreover, in the case of supers onic turbulence, the c l umps are expected to be 



by la r ge-scale compressiv e motions (IHunter fc Fleck 



2003 



Klessen et al 



1982 : 



Ballesteros-Paredes et al. 



ormea 



1999a . 



20051 ). suggesting that the turbulen ce within the clumys is similarly 



Vazquez- Semadeni et al 



(120081) 



likely to have a strong compressive component. Indeed, 
presented a numerical study of turbulent clouds driving showing that, even with continuous 
random driving at large scales, there is a trend for the clumps to contain a net convergence 
of the flow. This result contradicts the hypothesis that the non-thermal motions consist 
of isotropic, small-scale turbulence. After all, the clumps must be formed by a globally 
converging flow, as dictated by the continuity equation. 

n addition, numerical simulat i ons o f cloud formation in the diffuse atomic ISM 



star f ormation in turbulent, self-gravitating clouds (e.g . 



2000 



Heitsch. Mac Low &: Klessen 



2001 



Klessen 



2001 



(e.g. 


VazQuez-Semadeni et al. 


1995. 


1996: 


Passot et al. 


1995 




Ba 


lesteros-Parec 


es et al. 


1999a 




Kovama & Inutsuka 


2002: 


Audi 


; & Hennebelle 


2005 




Heitsch et al. 


2005 




2006. 


2008 




Vazquez-Semadeni et al. 


2ooe 




2007 




Hennebelle & Audit 


2007 


) and of 



Klessen. Heitsch fc Mac Low 



Bate. Bonnell. &: Bromm 



2003 



Vazquez-Semadeni et al 



Federrath et al. 



2008 



2005 

3~si 



Jappsen et al. 



2005 



Vazquez-Semadeni et al. 



2007 



20091 ) show that the velocity fields are in general organized at 



all scales, exhibiting a continuity from the large scales outside the clumps all the 
way to their interiors. Specific ally, the simulations of molecular cloud formation by 



Vazquez-Semadeni et al. 



(120071 . hereafter Paper I), showed that, as the cloud's mass is 
consumed by the formation of collapsed objects ("stars"), the cloud continues to accrete 
mass from its surrounding, more diffuse ("atomic") environment, the whole process 
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am ounting to a "mass cascade" , analo gous to a turbulent energ y casc ade, as proposed 



by 



Newman fc Wassermanl (119901 ) and 



Field. Blackman. fc Ketd (2008). Observationa. 



eviden ce supporting this scenario has been recently provided by 



Galvan-Madrid et al. 



too% . 



In this paper, we investigate the consequences of this hierarchical fragmentation 
scenario on the formation of low- and high-mass star forming regions, taking advantage 
of the variety of star forming regions appearing in the main simulation of Paper I. This 
study was motivated by a careful observation of the mechanism through which the clouds 
are assembled in that simulation. Specifically, we show that the region in the simulation 
with physical conditions resembling those of observed massive star-forming regions forms as 
the result of the large-scale collapse of a large complex, thus suggesting that such regions 
may be in generalized gravitational collapse rather than in a quasi - static state supported 



Hartmann &: Burkert 



(120071 ) for the Orion 



by turbulence, as proposed, for example, by 
molecular cloud. 

The plan of the paper is as follows. In §5] we briefly recall the main parameters and 
evolutionary features of the numerical simulation, and identify the regions we investigate. In 
^we measure the evolution of the total (dense gas + "stars") and "stellar" masses, and the 
instantaneous star formation rate (SFR) in both the isolated, peripheral, low-mass region 
(LMR) and in the central, global-collapse HMR. In ^ we then compare our results with 
available obervations. In §4.11 we compare the physical conditions and the star formation 
rates, showing that these indicators are reasonably similar to those of actual star forming 
regions of each type. Next, in §4.21 we perform a survey of the cores in our HMR, and 
compare the distributions o f their sizes, masse s and densities with those of a recent survey 



of the Cygnus X region by iMotte et al. 



(I2OO7I ). showing that the distributions are very 



similar. Finally, in ^ we discuss the limitations of our study and present a summary. 
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2. The numerical model 
2.1. Parameters 

In this paper, we analyze data from the simulation labeled L256At'0.17 in Paper I. We 
refer the reader to that paper for details. Here we just mention that it is a smoothed part icle 



hydro dynamics (SPH) simulation, performed with the GADGET code (ISpringel et al. 



20011 ). including self-gravity and parameterized heating and cooling that imply a thermally 
bis table medium. The s imulation uses 3.24 x 10^ SPH particles, and uses the prescription 



by 



Jappsen et al. 



(120051 ) for forming "sink particles" at densities above a certain density 
threshold, which we take as rithr = 3.2 x 10® cm~^. Once formed, the sinks can continue 
to accrete mass from their surroundings. Note, however, that the sinks are not sufficiently 
resolved in mass to be considered individual stars, and should be considered stellar clusters 
instead. The total mass in the simulation is 5.26 x IO^Mq, the size of the numerical box 
is 256 pc, and the initial conditions are uniform, with T = 5000 K and a mean density 
of n = 1 cm~^, in which two cylindrical flows are set up to collide at the center of the 
simulation. The flows have the same density and temperature as the rest of the box, but 
oppositely directed velocities of 7.536 km s~^, corresponding to a Mach number of 1.22. 
The flows have a length of 112 pc and a radius of 36 pc. 

It is worth noting that the cylindrical inflows are much narrower than the length of the 
numerical box, having a diameter of only 1/4 of the box length. This allows for interaction 
between the cloud and its environment, and also frees the cloud from boundary effects. 
Also, the length of the inflows is slightly shorter than the box length, so that the inflows 
are entirely contained within the box. This is because the boundaries are periodic, and thus 
the inflows cannot enter from the outside. This implies also that the total mass in the box 
is flxed. 
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2.2. Evolution 

The simulation we consider in this paper was designed to represent the formation 
and initial evolution of a dense, cold cloud out of the collision of streams of warm 



phase ( 


Hennebelle & Perault 


1999 




Wak 


er & 


i^'olini 


2000 


; I 


<ovama & Inutsuka 


2002 




Heitsch et al. 


2005; 


Audit & Hennebelle 


2005; 


Heitsch et al 


2006; 


VazQuez-Semadeni et al. 


2006. 


2007; 


Hennebelle & Audit 


2007 


), forming a large, flattened, and cold atomic cloud. 



which is turbulent because of a variety o 



warm components (iKoyama fc Inutsuka 



instabilities in the inter 



2002 



Heitsch et al 



2006 



ace between the cold and 



Vazquez-Semadeni et al. 



20061 ). The turbulence is moderately supersonic with respect to the cold gas. 



Because the cloud is colder and denser, it readily becomes self-gravitating and begins 
contracting and fragmenting, so that the large atomic cloud becomes a large molecular-like 
cloud complex. The densest fragments begin to undergo small-scale, local collapse events 



( iBanerjee et al 



20091 ). These isolated collapse events occur in a peripheral ring, which 



for ms because the collapse of a finite sheet-like cloud proceeds from the outside in, as shown 



by 



Burkert fc Hartmann 



(120041 ). Also, the local collapses occur before the bulk of the cloud 
completes its global, large-scale collapse, and the clumps and the stars they form are part 
of the global collapse of the complex. In the remainder of the paper, we show that the site 
at which the global collapse finally converges has physical conditions very similar to those 
of HMRs, while a typical local-collapse region has physical conditions resembling those of 
LMRs. 

The large atomic cloud formed by the flow collision, which later becomes a denser 
"molecular cloud''^ complex permeated by an atomic substrate, has a flattened shape 



^The simulation does not follow the chemistry, but we refer to gas with density and 
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during most of its evolution. The large-scale gravitational contraction begins at t ~ 10 
Myr, and star formation begins at t ~ 17 Myr. The initial radius of the peripheral ring is 
~ 20 pc. By t = 23.4 Myr, the global collapse is completed (the ring has shrunk to nearly 
zero radius), although the turbulence in the cloud has caused the infall motions to have a 
random component, so that the collapse center spans several parsecs across. Animations 
showing the large-scale evolution of this simulation can be found in the electronic edition of 
Paper I. 



It is important to note that in Paper I we applied a prescription by .Franco. Shore fc Tenorio-Tagle 



( 119941 ) to estimate a posteriori the time at which the total stellar mass formed in the 
simulation would imply the presence of enough massive stars to reionize the whole 
complex. We found this time to be roughly 3 Myr after the onset of star formation, 
or roughly at t ~ 20 Myr. However, simulations including ionization-heating feedback 
(Vazquez-Semadeni, Colin & Gomez, in prep.) suggest that, while the feedback is capable 
of disrupting the small-scale clumps, it does not prevent the large-scale collapse to form 
the massive region, since the material that falls into the latter comes fro m much larger 



distances than the lo c al spheres of influe nce of the stellar sources (see also 



Dale et al. 



2005 



Dale fc Bonnell 



2008 



Peters et al. 



20091 ). In addition, the cloud continues accreting mass 



from the warm medium even while it is actively forming stars. Thus, it is safe to assume 
that the central massive region would form even if we had included stellar feedback in the 
present simulation. 



temperature typical of molecular g clS dbS "molecular" . 
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2.3. The regions 

In this paper, we focus on two different regions of "star" (i.e., sink particle) formation. 
As a representative example of an LMR, we consider the first region to ever form stars in 
the peripheral ring starting at t = 17.3 Myr, which involves relatively small mass in both 
gas and sinks. Other LMRs in the simulation are similar to this one. As an example of an 
HMR, we consider the central region, where the global collapse converges at t = 23.4 Myr. 
We choose this time because it is the time of maximum compression of the material (see 
the animations corresponding to Fig. [2] in the electronic edition). This region involves a 
substantially larger amount of mass than the LMR, and is the only such region formed in 
the simulation, because the total mass in the complex is only ~ 5 x 10^ M© at the time 
we explore this region. Thus, this is our only choice for this type of region. Presumably, a 
more massive simulation could form a larger number of such regions. 

Figure [T] shows a column density map of the LMR at t = 19.1 Myr. The image is 
integrated over 16 pc in the x-direction. In the electronic version of the paper, this figure 
corresponds to an animation showing the evolution of this region from t = 16.6 to 19.9 
Myr. The dots show the newly formed sink particles. Figure [2] shows images of the central 
8 parsecs of the simulation (the HMR), at t = 24.2 Myr, viewed edge-on (left) and face-on 
(right). In the electronic version, this figure shows animations of the same region from 
t = 22.6 to 25.2 Myr for both views. In the images and animations of Fig. [2] we do not show 
the sink particles because a large number of previously formed particles is already present 
in the region, and it is not easy to see the particles formed there. 



Fig. 1. — Column density plot of Cloud 1 in the y-z plane at t = 19.1 Myr, integrating over 
the central 16 pc along the x direction. The dots show the stellar objects (sink particles). 
The electronic version of this figure shows an animation of this region from t = 16.6 to 19.9 
Myr. The column density is in code units, which correspond to 9.85 x 10^^ cm~^. 
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3. Physical conditions of the "low"- and "high"-mass regions 

We begin our comparison study between the two regions by measuring the evolution 
of physical and star-forming properties of the "clouds" in these regions. To this end, we 
interpolate the SPH data of the regions shown in Figs. [T] and [2] onto a 256^ grid in order to 
manipulate the data with the IDL software. We then define the clouds as connected sets 
of pixels with densities above nth = 500 cm~'^. We refer to the resulting clouds as "Cloud 
1" (CI) in the LMR, and "the Central Cloud" (CC) in the HMR. In Fig. E] we compare 
the evolution of the physical properties of the two clouds over a few Myr. Note that the 
starting times (to) for the evolutionary plots are different in the two cases, being to = 17.27 
Myr for CI and to = 22.58 Myr for the CC. The thick lines correspond to the CC, and the 
thin lines correspond to CI. The left column shows the evolution of the density, size and 
velocity dispersion for the two clouds, where the size is calculated as r = 
being the volume of the cloud. 

From these figures we see that the CC clearly has larger density, size, velocity 
dispersion, and mass than CI, at least during the time spans shown in the figures (CI later 
acquires more mass and becomes more similar to the CC). Although the density of the CC is 
only a factor of 2 larger than that of CI on average, we see that its size ranges between half 
and one order of magnitude larger than that of CI, and similarly for i ts velocity dispersion , 
which is on the order of a few km s~^, as is the case for HMRs (e.g., iBeuther et al.l 120071 ). 
Similarly, the CC's mass ranges in the thousands of solar masses, while that of CI initially 
ranges in the hundreds, although, as mentioned above, this region is growing in time, and 
approaching the conditions of the CC at late times. Note that the gas mass in the CC is 
initially less than the sink mass because we include the mass of pre-existing sinks in the 
region, and not only the sinks formed by the CC. We do this because we next compute with 
the virial mass of the system, which must be compared with the total gravitational mass. 
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including the gas and all sinks in the region. 



We compute the virial mass of each region applying the standard formula 

M..^210(ii](.i^]'Me. 



pc J ykm s~^ 



(1) 



[see, e.g., 



Caselli et al. 



2002 



Tachihara et al. 



2002 



Klessen et al. 



20051 ). It is noteworthy 



that the CC, even though it is in an extremely dynamic infalling state, has a ratio of 
mass to virial mass M/Mj of almost exactly unity during most of its evolution. This is 
precisely what is expected for an object undergoing gra vitational contraction, and has 



been observed previously in n umerical simulations (e.g 



Klessen &: Burkert 



2000 



2001 



Vazquez-Semadeni et al. 



20071 ). Instead, this ratio is completely arbitrary in the case of 



clouds or clumps formed purely by turbulent compressions. 

Cloud 1, on the other hand, starts with a ratio of nearly 5, but it steadily approaches 



unity as it accelerates in its contraction. Note th at such values o 



uncommon in observed low- mass clouds (see, e.g.. 



Tachihara et al. 



■ this ratio are not 



2002 



Morata et al. 



20051 ). Presumably, the relatively large initial values of the ratio indicate that at its initial 
states, the cloud is only beginning to decouple from the global flow, and starting its own 
collapse. 



We furthermore measure the sink formation rate (SiFR) in both clouds, defined as 



SiFR 



At 



(2) 



where AMgink is the increment in sink mass in the region during the time interval At = 0.13 
Myr between successive data dumps. Analogously, we define the "specific SiFR" {sifr) as 

SiFR 



sifr 



(3) 



M ' 

^ " gas 

where Mg^s is the instantaneous dense gas mass in the region, again defined as the mass of 
gas with densities n > 500 cm~'^. Note that, for simplicity, we perform these measurements 
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directly on the SPH data, rather than on their interpolation on a grid. Note also that we 
use the increment in the mass of the sinks rather than in the number of sinks because, 
as mentioned in §21 our sinks are not sufficiently resolved in mass, and furthermore they 
continue to accrete mass after they form, as shown in Fig. |H This figure shows the mass 
evolution of the most massive sinks formed in the simulation up to t = 20 Myr {left panel), 
the most massive sinks formed in the central 8-pc region for 23.9 < t < 25 Myr {middle 
panel), and those formed outside of the central 8-pc region during the same time interval 
{right panel). 

Both the SiFR and the sifr are shown in Fig. O It is clearly seen that again the SiFR 
of the CC is roughly one order of magnitude larger than that of CI. However, it is seen 
that the specific star formation rate is similar for both clouds, suggesting that the larger 
SiFR of the CC is due exclusively to the larger amount of mass available for collapse in this 
cloud. In any case, if one assumes a universal IMF, then the larger SiFR of the CC implies 
that it has a ten-fold larger likelihood of forming massive stars, in correspondence with its 
more massive and violent status. 



4. Comparison with observations 

In this section we compare the mean physical properties of the massive clump in the 
CC with those of a typical massive-star forming clump in the Orion Molecular Cloud, 
namely the well-known clump OMC-1. Next we compare the statistical distributions of 
masses, sizes and densities of the dense cores in the Cygnus-X North region, reported by 



MotteetaL 



(120071 ) with the corresponding structures within the CC. 



It is important to note that, since we are making these comparisons against data from 
two different molecular clouds, we are making the implicit assumption that the massive-star 
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forming clumps in GMCs are similar in different clouds, constituting a general class of 
objects. This is, however, a reasonable assumption, since these clumps are often discussed 
in general terms, re gardless of the GMC to which they belong (see, e.g., the discussion in 



Beuther et al. 



2007, sec. 2.1). 



4.1. Physical properties and star formation rates 



We begin by considering a typical HMR, namely OMC-1, the gas behind the 
Orion Nebu l a Clu ster (ONC) in the Orion Molecular Cloud (OMC). According to 



Bally et al. 



(119871 ). the total gas mass in the clump known as OMC-1 is ~ 2200 Mg 







and the typical observed (i.e., ID) rms velocity is ~ 2 km s ^, with an angular size 



~ vT08 arcmin (se e the first entry in their Table 1). Assum ing a distance of 



400-500 pc to OMC-1 flGenzel et al 



1981 



Sandstrom et al 



20071 ) ■ this value of A9 



implies a physical size ~ 1.2-1.5 pc for this cloud. Thus, the mean density in OMC-1 is 
~ 22OOM0/(1.35 pc)^ ~ 900 Mopc'^ ~ 1.54 x lO'' cm~^ where for the last term we have 
assumed a mean particle mass of 2.36. 

The above numbers can be compared with the physical conditions in the CC in our 
simulation. To do this, we need to select the gas within this region with similar mean 
density. We have found that this is accomplished by using a density threshold of 4500 cm~'^ 
to define the region, which gives us a mean density of precisely 1.54 x 10^ cm~^. 

Next, we apply a clump-finding algorithm based on identifying connected sets of 
grid points whose densities are above this threshold, allowing us to identify the most 
massive clump within the CC. This region c an be considered th e equivalent to what is 



often referred to as a "cluster-forming core" (iLada fc Lada 



20031 ). However, we refer to it 



as a "clump". Once the clump is defined, we can straightforwardly measure its physical 
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properties. We find a gas mass Mgas ~ 3015 M©, a size r ~ 1.14 pc (computed through 
the formula r = {SV/AnY^^, where V is the volume of the clump), and a three-dimensional 
velocity dispersion a ~ 4.74 km s~^, implying a one-dimensional velocity dispersion of 
~ 2.74 km s~^. We thus find that there is very good agreement between the physical 
conditions in OMC-1 and those of an equivalent region within the CC in our simulation. 

Moreover, from Fig. [5] we see that the CC is characterized by a SiFR ~ 10~^Mq yr~^ 
on average. Noting that our density threshold for the formation of a sink, of 3.2 x 10^ cm~^, 
is comparable to the density of a dense core, we should consider that the star formation 
efficiency (SEE) within our sinks is less than 100%. Assuming an SEE ~ 30- 50%, we 
estima te a star formation rate (SER) ~ 3-5 x IO'^Mq yr^^. On the other hand 



Tobin et al. 



( 120091 ) report on 1613 stars in the ONC. Taking this number as a proxy for the total stel 



prod uction of this region, and a mean stellar mass of 0.3 Mq (iHillenbrand fc Carpenter 



ar 



20001 ). this ir nplies a tota . 



is < 2 Myr (IHillenbrand 



stel 



ar mass of ~ 500Mq. The estimated age spread of the cluster 
1997f ). We thus infer SER > 2.5 x lO'^M© jt~\ in very good 
agreement with that of our CC. 

With these data, we can also calculate the SEE per free-fall time, SEEgi Specifically, 
writing 

where Mstars(^ff) ~ SER x tg, with SER being the instantaneous star formation rate, 
estimated above, and tg the free-fall time. At the mean density of 1.54 x 10'^ cm~^, and a 



^This quantity is frequently referred to as SER per free-fall time in the literature. How- 
ever, it is actually the specific SER, integrated over a free-fall time (e.g. iKrumholz &: McKee 
2OO5I : Krumholz fc Tanll2007l ). 
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temperature of 10 K, we have tg ^ 0.27 Myr. We thus find 

135M0 



SFEfl 



0.043. 



315OM0 

Interestingly, this value is fully consistent with that obtained by 



(5) 



Krumholz fc TanI (120071 ) 



(SFEff ~ 0.03-0.09) for the ONC. We discuss this resuh further in CT 

Finally, we can also perform a compa rison between the CC and OMC-1 at the level of 



the core formation rate. 



IkedaetaL 



(120071 ) report a total mass of dense cores in the OMC-1 
regio n of ^ 800 Mf^. Assu ming that these have typical lifetimes between 2 and 5 x 10^ yr 



[e.g. 



Hatchell et al. 



120071 ). we find a core formation rate (CFR) between 1.6 and 4 x 10 M© 
Myr~^, which again compares very well to the mean value of the SiFR mentioned above. 
Note that the estimated CFR for OMC-1 is larger than the estimated SFR in the ONC by 
a factor of a few, a reasonable result if the efficiency of star formation within the cores is 
smaller than unity, and consistent with the factor we have used above for estimating the 
SFR from the SiRF. 

We conclude that the physical conditions and SFR in a well-known massive-star 
forming clump agree to within less than factors of 2 with the corresponding features in the 
CC of our simulation. 



4.2. Statistics of cores in the Central Cloud 

As a further test of the similarity of the CC with massive-star forming clumps, in this 
section we c o mpar e the statistics of the dense cores within the CC with those reported by 



Motte et al. 



(120071 ) for the Cygnus X molecular complex. These authors have performed a 
1.2 mm continuum unbiased survey of the massive-star forming clumps within the Cygnus-X 
complex, and reported the masses, sizes and densities of the cores within the clumps, 
and of the clumps themselves. Specifically, they compiled a large dataset of 129 massive 
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dense cores, of sizes ~ 0.1 pc, masses Mi. 2 mm = 4-950 Mq, and densities n ~ 10^ cm~^, 
listed in Table 1 of that paper. These data are very well suited for comparison with the 
corresponding structures in the CC of our simulation. In order to perform the comparison, 
we identify dense cores within the gridded version of the CC in our simulation (cf. Fig. 
[2]) at time t = 23.9 Myr, at which the global collpase of the entire cloud complex in the 
simulation is being completed (i.e., the clouds that had initially formed at the peripheral 
ring in the simulation reach the center at this time). 



To define the cores, we apply the same clump-finding algorithm used in §4.11 but in this 
case we use it iteratively and with higher threshold densities to define an ensemble of cores. 
Specifically, we consider a set of four thresholds, at n = 10^, 3 x 10^, 10^, and 3 x 10^ cm^^, 
and include all of the resulting data points into a single dataset. We do this in order to have 
a sufficiently large sample, and to allow for a wider dynamic range in the physical properties 
of the cores, since our previous experience shows that objects selected with a given threshold 



have mean den sities that differ by less than an order of magnitude 



threshold (e.g., 



Vazquez-Semadeni. Ballesteros-Paredes fc Rodriguez 



proce dure differs from other algorithms such as CLUMPFIND (jWilliams. de Geus fc Blitz 



rom the value of the 



199711. Note that our 



19941 ) in that, in our case, clumps defined at a lower threshold may contain several clumps 
defined at a higher threshold, while CLUMPFIND, for example, does not identify such 
"parent" structures as a clump. We choose this procedure for simplicity, and because larger 
clumps defined at l ower thresho l ds ma y be bound structures themselves that should be 
taken into account. 



Motte et al 



(|2007f) themselves used a combination of source extraction 



analysis and GAUSSCLUMPS (Stutzki fc Glisten 



1990). 



This procedure leaves us with 39 cores, 38 of which have M > 4Mq. For each 
core, we measure its mass, mean density, and velocity dispersion, and estimate its size 
as r ^ (3V"/47r)^'''^, where V is its volume. The distributions of these properties can be 
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compared with those reported by iMotte et aLl (120071 ) for any of the regions within the 
Cygnus-X complex. As a representative example, we choose the data for the Cygnus-X 
North region, which contains 72 cores (see their Table 1). We show the distributions of 
core properties for the Cygnus-X North region in the histograms presented in Fig. [6] by 
the dotted lines. Superimposed on these histograms, the solid lines show the corresponding 
distributions for the 38 cores more massive than AMq in the CC of our simulation. We 
see that the distributions match extremely well for size , mass, and i nean density. No 



comparison is made for the velocity dispersions because 



Motte et al 



(120071 ) do not report 



them, but in our sample, they range from ~ 0.2 to ~ 3 km s . 

A final comparison that can be made is at the "clump" level. This is analogous to 
the comparison we made with OMC-1, but comparing against the mean properties of the 



clumps identified by 



Motte et al. 



(120071 ). They report the typical sizes, masses, and densities 
of all clumps to be r ~ 0.68 pc, M ~ lOOOM©, and n ~ 1.4 x 10^ cm~^ respectively (see 
column 3 in their Table 4). Again, this is comparable to the clump we obtain thresholding 
our CC data at nth = 4500 cm~'^, for which we find a mean density of n 1.54 x 10^ cm~^, 
a size r ^ 1.14 pc, and a mass of 3015 Mf^ (c f. ^4.11) . We see that our clump is somewhat 



more massive and extended than 
of 2-3 of the latter. 



Motte et al. 



(j2007l )'s mean clump, but still within factors 



5. Discussion and caveats 

5.1. The star formation efficiency per free-fall time 

The results from §4.11 have a number of interesting implications. We note that the 
SF rates in both the CC in our simulation and in OMC-1 in Orion are relatively large, 
~ 3 X 1O~^M0 yr^^. This can be translated into a specific SFR (or inverse SF, or gas 
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depletion, timescale, tg^), given by 



SFR 



(6) 



where Mgas is the gas mass in the cloud and M* is the stellar mass. We obtain tgr ~ 10 
Myr for both the CC and OMC-1. T his differs strongly from the much longer timescales 



tsF ~ 2 Gyr for extragalactic GM Cs (IBigiel et al. 



GMCs flMcKee fc Wilhams 



20081 ) ■ and tgF ~ 300 Myr for Galactic 



19971 ). Said another way, if one were to extrapolate (linearly 



with mass) the SFR of these clumps (the CC or OMC-1) to the entire molecular gas 
mass of the Galaxy of IO^Mq, one would find a total Galactic SFR of ~ 3 x IO^'^Mq 
yr^^ X IO^Mq/SOOOMq = 10 M(d yr~^, or roughly 50-10 times larger than the observed 



SFR. This is the well known 



Zuckerman fc Palmen (119741 ) conundrum. However, also in 



section §4.H w e showed that the SFEff of both regions is small, and fully consistent with 
the estimate of iKrumholz fc TanI (20071 ) for O NC, and the general result of a low SFEff for 



the cores in the survey by 



Evans et al. 



(120091). 



These results then raise two important questions. The first is how to reconcile the 
much larger specific SFR of regions like the CC or OMC-1 with that of entire molecular 
clouds, or even the total molecular gas content of the Galaxy. The large disparity in 
tgp between individual cluster-forming clumps and large-scale objects like GMCs can be 
understood if star formation is a spatially intermittent phenomenon. That is, it does not 
occur everywhere in the molecular gas, but only at special locations, characterized by the 
highest densities. Thus, the small fraction of the mass that is forming stars must have a 
much larger specific SFR than the global average. This is analogous to the intermittency of 
the energy dissipation rate in incompressible turbulent flows, which is known to occur only 
at sc attered locat ions, where the velocity gradient is maximized by vortex stretching (see. 



e.g. 



Frisch 



19951) 



The second question is: how can such an actively star forming region have such a small 
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value of the SFEg. The theory by iKrumholz fc McKed (120051 ) originally derived such low 
values for the SFEg for clumps in near virial e quilibrium supported by turbulence (similarly 



McKee fc Tan 



20031 ). However, the CC in 



to the underlying assumption of the model by 
our simulation is far from being in equilibrium, as shown in the animations corresponding 
to Fig. 

The origin of the small SFEfj value for the CC in our simulation can be understood by 
rewriting eq. (jl]) as 

SFEff = §^ tff, (7) 

where Mtot is the total mass (gas + stars) in the cloud. From this expression we see that a 
low value of the SFEg can be obtained through a large cloud mass and/or a high density 
(via tff), besides the obvious possibility of a small SFR. This appears to be the case of the 
CC, which is rapidly increasing its mass through accretion of the surrounding gas at a rate 
that apparently overwhelms its consumption rate by SF. Indeed, Fig. [7] shows the evolution 
of the total gas mass in the central 8-pc box of the simulation in which the CC is contained. 
The gas mass is seen to increase from ~ 4OOOM0 to ~ I2OOOM0 in 3 Myr, corresponding 
to a mean accretion rate onto the region of ~ 2.7 x 10^'^ Mq yr^^, or roughly an order of 
magnitude larger than the SFR in the CC. 

Moreover, because of the ram pressure (or, equivalently, weight) of the accreting 
material, the density is also larger than what would be expected if the clump were isolated 
and in hydrostatic equilibrium. That is, the actual collapse timescale of a region such 
as our CC and its surroundings is longer than the free-fall time of the central dense core 
alone, because the mean density of the combined system is lower than that of the core, and 
the material is continuously being replenished at the central regions. Thus, if the infalling 
material is unseen for any reason (e.g., thresholding, background substraction, termination 
of the tracer excitation, etc.), then the central core will appear to have a shorter free-fall 
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time than the actual timescale of the collapse. This can be incorrectly interpreted as an 
isolated, equilibrium core with a lifetime much longer than its local free-fall time and thus 
with an SFR much slower than that given by the free-fall rate. 

Finally, it is worth noting that this scenario the n naturally explains the large final 



Lada fc Lada 



20031 ). without having 



efficiencies of cluster-forming clumps (~ 30-50%; e.g., 
to invoke lifetimes for these clumps that are much larger than their free-fall times. Also, 
in this accreting scenario, we see that a low local value of the SFE^ does not guarantee a 
globally low SFR. Our results show that a region can have a large SFR and still have a 
low SFEfj. The attainment of a globally low SFR seems to require a highly scattered and 
intermittent occurrence of the star-forming regions. 



5.2. Neglect of stellar feedback 



The simulation we have analyzed in this paper neglects two important processes 
existing in real molecular clouds, namely stellar feedback and magnetic fields. Concerning 
the former, it is generally believed that it is a key mechanism in keeping a low final SFE 
in molecular clouds and their clumps. However, at present it is still unclear whether this 
self-regulation of SF occurs abruptly, after a brief period of a high SFR, dispersing the 



star- 



'orming region and quenching SF a . 



Whitworth 



1979 



Elme green 



1983 



together locally 



Elmegreen 



( the so-called "rapid SF" scenar io; 



2000 



Elmegreen 



2007 



Larson 



1987 



Franco. Shore fc Tenorio-Tagle 



1994 



Hartmann. Ballesteros-Paredes. &: Bergin 



200l|), or 



whether it occurs continuously over longer times, maintaining the region in a near-equilibrium 



state, allowing it to last for seve ral free-fall times, w. 



so-cal 



2006 



ed "slow SF" scenario ; e.g. 



Li fc Nakamura 



2006 : 



lile n i aintaining a low SFR (the 



Krumholz fc McKee 



Krumholz fc Tan 



20071) 



2005 



Krumholz. Matzner fc McKee 
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The behavior of our simulation is consistent with the "rapid SF" scenario, in which a 
given region forms stars at a high rate over a short period of time and then their feedback 
disrupts the cloud, possibly assembling a new star forming region at a new location , through 



the "collect and collapse" mechanism (lElmegreen fc Lada 



1977 



Elmegreen 



20071 ). Since 



the large complex and its fragments in our simulation begin their gravitational contraction 
several Myr before the first star-forming collapse events occur (cf. §2.2p . it is natural to 
expect that the initial stages of SF in each region will occur in a gravitationally contracting 
environment, and to be characterized by a large SFR, until enough stars hav e been form ed 
that they can strongly affect, and possibly disrupt their parent clump (e.g.. 



1983 



Larson 



1987 



Franco. Shore fc Tenorio-Tagle 



199J). 



Elmegreen 



Under this scenario, our simulation and the star-forming regions formed within it 
represent the assembly and early star-forming epochs of a region, although their long-time 
star formation efficiencies are too large (~ 60-70% after ~ 10 Myr; cf. Paper I) because they 
lack the subsequent disrupting effect of the stellar feedback on the clump in which the stars 
form, and so SF continues unimpeded in each region until the gas is exhausted. However, 
this is a shortcoming that is only expected to become important later into the evolution 
of each region, when the accumulated number of stars becomes large enough to exert a 
str ong influence on the cloud. As mentio ned in §2.21 in Paper I we applied a prescription 



by 



Franco. Shore fc Tenorio-Tagle 



(119941 ) to estimate this time, which we foun d to be ~ 3 



Myr after the onset q : 



2005 



Dale fc Bonnell 



SF, although new simu lations including stellar feedback (jPale et al. 



2008 



Peters et al. 



20091 . Vazquez-Semadeni et al., in prep.) suggest 



that the effect of feedback is more complicated than that, and should be applied clump by 
clump. But, in any case, our simulation should be accurate before the time of disruption 
by feedback, for each new star-forming region formed. 



On the other hand, the evolution of our simulation is inconsistent with the "slow" 
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scenario of star formation (e.g.. iKrumholz fc McKedl2005l : iKrumholz fc Tarul2007l ). in which 
clouds are assumed to be maintai ned in qu asi-st atic e q uihbrium by the turbulence inj ected 



by st el 
20061 ): 



ar sources, as m odelec 



Li fc Nakamura 



(e.g.; 



Matznerl (e.g., 



2002); 



Krumholz. Matzner &: McKed (e.g., 



20061 ). Our simulation, although incapable of following the 



stages in which SF interacts strongly with the infall process onto the clouds, suggests that 
they should be already contracting in general by the time they begin forming stars, rather 
than being in a quasi-hydrostatic state. Besides, the studies suggesting such a hydrostatic 
state have been performed either in closed numerical boxes or with simplified geometries, 
so that the disruption of the cloud by its stars has been made overly difficult. Moreover, it 
appears unlikely that the effect of the stellar energy injection, which occurs at small scales, 
can somehow organize itself to pr oduce the large-scale, dipolar character of the principal 



velocity component in the clouds (IHeyer fc Brunt 



20071 ). Simulations in which the cloud is 



embedded in its diffuse medium and in full three-dimensional geometry, such as the one 
utilized here, but including the effects of stellar feedback, are necessary to properly address 
the issue (Vazquez-Semadeni et al., in prep.) 

Our results are consistent with recent suggestions, ba sed on comparisons betwee n 



simulations and observations, that molecular c . 



and clumps (e.g. 



Peretto. Hennebelle fc Andre 



ouds (e.g.. 



Hartmann fc Burkert 



20071 ) 



20071 ) may be in a state of gravitational 



collapse. If confirrned, the se suggestions point towards a return to the original suggestion by 



Goldreich fc Kwan 



(119741 ) that the observed linewidths in mole cular clouds are due primarily 



Zuckerman fc Palmed (119741 ) 



to gravitational contraction. This suggestion was dismissed by 
through the argument that this would imply a much larger average star formation rate in 
the Galaxy than observed. However, in the scenario of rapid SF the instantaneously high 
rates of SF are not a problem, if the star-forming clumps are soon disrupted by the stellar 
feedback before their entire mass is consumed by SF. 
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5.3. Neglect of magnetic fields 



Our simulation also neglects the influence of magnetic flelds, and as such it should be 
considered as representative of the evolution of the supercritical parts of the clouds, which 



are the ones that can undergo global collapse. For t 
but nonzero magnetic flelds can retard the collapse ( 



lese regions, the presence of weak 



Ostriker. Stone fc Gammie 



2001 



Heitsch. Mac Low fc Klessen 



Vazquez-Semadeni et al. 



2001 



20051 ). In any case, if other parts 



of the clouds are sub critical and are thus held up against their self- gravity by the magnetic 
fleld, th en they must con tribute only a small (or nearly null) fraction of the SF in the 



Galaxy (Elmegreen 



( ILada fc Lada 



20071 ). since it is known that most stars form in cluster-forming regions 



20031 ). while a subcritical region is expected to be a site of isolated, low-mass 
star formation. Of course, if a large fraction of the cloud's mass is supported by the fleld, 
then the global SFE for the cloud will certainly be decreased, but the supercritical regions 
will still be responsible for the majority of the SF activity of the cloud, and those parts 
should have similar SFEs as our simulation. Moreover, as discussed in §5.11 it is precisely 
this kind of structure that is necessary to have a large fraction of the molecula r mass not 
parti cipating in SF, while the remaining small fraction (a few percent; see, e.g. 



Kirk et al. 



20061 ) can be essentially free-falling, and accounting for most of the SF in the GMC. 



5.4. Clump properties 

Although the CC has been shown to reproduce several features of massive-star forming 
clumps, it does also exhibit certain differences. Most notably, the total mass in the CC 
contains roughly twice the mass in sinks than in gass at t — to = (cf. Fig. [3l upper right 
panel). This is contrary to the situation in the ONC-OMC-1 complex, where the gas mass 
is roughly 4 times larger than the stellar mass. The excess sinks in the CC are "stragglers" , 
which were formed earlier in the peripheral ring, but that have fallen into the large potential 
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well of the whole complex together with the gas. This suggests that the global collapse 
responsible for the formation of the CC is somewhat too focused, perhaps as a consequence 
of the extremely smooth initial conditions used in the simulatio n. Indeed, experimentation 



20091) 



with more strongly fluctuating initial conditions in the velocity (iRosas-Guevara et al. 
shows that in this case the collapse is less focused, possibly causing the population of 
infalling sinks in the collapse center to be less numerous. Nevertheless, the presence of such 
i nfalling stars may explain the age spreads observed in regions of massive star formation 



flPalla fc Stahler 



1999 



Faustini et al 



20091) 



Note also that, due to the presence of the infalling sinks, it is possible for both the 
gas and sink mass in the upper right panel of Fig. [3] to decrease simultaneously. This is a 
consequence of the gas being consumed to form new sinks, but with the SiFR not being 
able to compensate for the departure of the infalling sinks which, contrary to the situation 
of the gas, simply pass through the region, while the gas shocks and stagnates there. 

Finally, it is necessary to remark that the dense cores necessarily lie at the limit 
of the scales resolved by the simulation. Thus, the velocity dispersion within them is 
most certainly reduced by numerical diffusion. However, the fact that the gravitational 
contraction is hierarchical, and starts at the largest scales of our cloud, implies that the 
collapsing nature of the clumps cannot be an artifact of the numerical diffusion. 



6. Summary and Conclusions 

In this paper we have presented numerical evidence that the physical conditions in 
low- and high-mass star forming regions ( "clouds" ) can arise from hierarchical gravitational 
collapse. The former regions arise from small-scale fragments in the collapse (which however 
occur first, because they originate from local high-amplitude density fluctuations which 
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have a shorter free-fall time), while the latter may appear when the global, large-scale 
collpase is completed. The local collapse events do not exhaust the gas in their regions 
because a) the consumption time is relatively long, ~ 5 Myr (cf. Fig. right panel), and 
b) because the clouds continue to accrete mass from their atomic surroundings throughout 
their evolution. Thus, the isolated low-mass clouds eventually collide to form a high-mass 
complex. In this scenario, velocity dispersions are caused primarily by infall motions rather 
than random turbulence, but its hierarchical nature may explain why massive cores tend 
to have larger velocity d ispersions than low-mass ones at the same size (see, e.g.. Fig. 10 



of 



Garay fc Lizano 



19991 ). since the former only begin contracting when they decouple from 



the global flow, while the latter have been accelerating in their contraction for a longer time. 

The evidence we presented consisted of two parts. First, we analyzed the evolution of 
the physical properties (mass, mean density, size, velocity dispersion, and mass to virial 
mass ratio), as well as of the SFRs of two examples of the class of objects, which we called 
"Cloud 1" (CI) and "the Central Cloud" (CC). We noted that the latter, which forms as a 
consequence of the large-scale collapse of the cloud complex in the simulation, contains a 
clump with properties consistent with those of cluster-forming cores within GMCs, such as 
a gas mass Mgas ~ SOOOMq, a velocity dispersion of a few km s~^, and a star formation 
rate SFR ~ 3 x lO~^M0yr~^, while CI is more reminiscent of low- or intermediate- mass 
regions, with Mgas ~ 400Mq, ~ 0.7 km s^^ (during its initial stages, as it later becomes 
more massive), and SFR ~ 3 x lO~^M0yr~^. 

Secondly, we performed a survey of the dense cores in the massive cloud (the central 8 
parsecs of the simulation at the time the large-scale collapse converges there), and compared 
their statist ical properties with those reported for the dense cores in the Cygnus X cloud 



complex by iMotte et al. 



(120071 ) ■ finding that the distributions of both samples match each 



other very well. At the slightly larger-scale, "clump" level, we also found good agreement 
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between the mean properties reported by those authors and the "clump" in that region. 
These results suggest that indeed the convergence of the global collapse is a good model of 
HMRs. 

This scenario, however, has the implication that the mechanism responsible for the low 
global Galactic SFR may be the result of SF being a spatially intermittent process, so that 
most of the molecular gas mass is not forming stars at any given time in the Galaxy, while 
those regions that are forming stars do so at a very high rate — a rate that in fact accounts 
for most of the SF activity in GMCs and the Galaxy. This, in turn, implies that the global 



specific SFR of the molec ular gas in the Ga . 



Krumholz fc Tan 



axy d oes not extrapolate linearly with mass (as 



20071 ) to estimate the specific SFR of a single 



assumed by, for example, 
SF region. 

Finally, we found, somewhat unexpectedly, that the massive-SF-like region in our 
simulation exhibits a low SFEg (~ 4%), in spite of having a large SFR of ~ 3 x 10~^Mq 
yr~^, and it being the site of the large-scale collapse of the cloud formed in the simulation, 
rather than a hydrostatic structure supported by turbulent pressure. We attribute this to 
the high density and large accretion rate occuring in this region as a consequence of the 
culmination of the large-scale collapse. In particular, the region accretes mass at an even 
higher rate than that at which it is forming stars, even though the latter rate is large. 
This physical scenario calls for a new class of theoretical models that take into account the 
sources and sinks for the mass of a molecular cloud. 
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Fig. 2. — Two views of the central 8-pc cubic region. Left panel: Column density integrated 
along the y direction. Right panel: Column density integrated along the z direction. The 
electronic version of these images show the evolution of these regions from t = 22.6 to 25.2 
Myr. 
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Fig. 3. — Evolution of the physical properties of the dense gas (n > nth = 500 cm~^) in 
the two star-forming regions anlyzed in the simulation. The thick lines refer to the "Central 
Cloud" (CC), and the thin lines refer to "Cloud 1" (CI). Left column: Evolution of the 
mean density {top panel), size {middle panel), and velocity dispersion {bottom panel). Right 
panel: Evolution of the dense gas mass and sink mass {top panel); evolution of the total 
(dense gas+sinks) mass and the virial mass {middle panel); evolution of the ratio of total 
to virial mass {bottom panel). Note that the starting times to for the evolutionary plots are 
different for the two regions, being to = 17.27 Myr for CI and to = 22.58 Myr for the CC. 
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sinks formed before t = 20.0 Myr inner (8pc)"^ 23.9 < t < 25.0 Myr not inner (8pc)^ 




Fig. 4. — Evolution of the most massive sinks formed in the simulation up to t = 20 Myr 
{left panel), the most massive sinks formed in the central 8-pc region for 23.9 < t < 25 Myr 
{middle panel), and those formed outside of the central 8-pc region during the same time 
interval {right panel). 




Fig. 5. — Evolution of the SiFR and the specific SiFR (si/r) for Cloud 1 (CI, dotted lines) 
and the Central Cloud (CC, solid lines). 
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Fig. 6. — Histograms of the size {top left), mass {top right), mean density {bottom left) 
and three-dimensional velocity dispersion {bottom right) of the cores within the Central 
Cloud of the simulation {solid lines) at t = 23.9 Myr, at which time the region acquires 
its most compressed state. The histograms are normal ized to their peak values. Dotted 



Motte et al 



(120071 ). except for the 



lines: same for the core sample in Cygnus X reported by 
velocity dispersion, since their observations were performed in the 1.2 mm continuum. The 
distributions are seen to match extremely well for size, mass and mean density. 
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Fig. 7. — Evolution of the total gas mass in the central 8-pc box of the simulation. The 
total mass in this region is seen to increase from ~ 4 x 10^ to ~ 1.2 x 10^ in 3 Myr, 
corresponding to a mean accretion rate of 2.7 x 10~^Mq yr~^. 
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